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Abstract 

Specific forms for both the Gibb’s and the complementary dissipation potentials were chosen such 
that a complete potential based multiaxial, isothermal, viscoplastic model was obtained. This model 
in general possesses three internal state variables (two scalars associated with dislocation density 
and one tensor associated with dislocation motion) both thermal and dynamic recovery mechanisms, 
and nonlinear kinematic hardening. This general model, although possessing associated flow and 
evolutionary laws, is shown to emulate three distinct classes of theories found in the literature, by 
modification of the driving threshold function F. A parametric study was performed on a specialized 
non dimensional multiaxial form containing only a single tensorial internal state variable (i.e., internal 
stress). The study was conducted with the idea of examining the impact of including a strain- 
induced recovery mechanism and the compliance operator , derived from the Gibb’s potential, on 
the uniaxial and multiaxial response. One important finding was that inclusion of strain-induced 
recovery provided the needed flexibility in modeling stress-strain and creep response of metals at low 
homologous temperatures, without adversely affecting the high temperature response. Furthermore, 
for nonproportional loading paths, the inclusion of the compliance operator had a significant influence 
on the multiaxial response, but had no influence on either uniaxial or proportional load histories. 
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1 Nomenclature 


Invariants 

n 

complementary dissipation potential 

4 

Gibb’s complementary potential 

F 

Bingham- Prager threshold function 

G,4 

normalized second invariant functions 

J 2 

second invariant of effective deviatoric stress tensor 

h 

second invariant of internal deviatoric stress tensor 

Hi, Hi 

invariant material functions 

Stresses 


Cauchy stress tensor 

Si. i 

deviatoric stress tensor 

E.-i 

effective deviatoric stress tensor 

°7 

internal state variables (stress-like) 

“•j 

internal (or back) stress tensor 

a ij 

deviatoric internal stress tensor 

K,K C 

drag stress and initial drag stress 

Y 

yield stress 

Strains 

a ■ • #1. f/? f T! 

> c ij > S; > Sj 

total, inelastic, reversable, and thermal strain tensors, resepectively 

Ay 

conjugate internal state variables (displacment-like) 

-4»> 

internal strain tensor 

/C 

drag strain 


yield strain 

Material Parameters 

Cijki 

elastic moduli tensor 

H,R 

hardening and thermal recovery material parameters 

fi,B 

material parameters 

0 

denotes the extent of strain induced recovery 

n,m,p,q 

material exponents 

g(G),f(F),x( K ),z(Y) 

material functions 

Efaj )< H(a y ), Z(T) 

material functions 

n 

coefficient of thermal expansion 

Miscellaneous 

T,T m ,T melt ,T 0 

applied, homologous, melting, and reference temperature, respectively. 

t 

real time 

T 

nondimensional time 

n,B 

nondimensional groups representing classes of materials 

gi,Q2,Q3 

distinct compliance operators 

6y 

Kroncker delta function 


indicates variable normalized by initial drag stress, k 0 

0 

Macauley Bracket 

Hv[ ] 

Heaviside unit function 

(° ) 

nondimensional time derivative notation 

(• ) 

time derivative (or rate ) notation 
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2 Introduction 

In a previous publication [1], we reported on the general framework of viscoplastic representations based 
on the thermodynamics of irreversible processes and the internal variable formalism (for isotropic metals 
in the small- deformation regime under isothermal conditions). Emphasis in reference 1 was placed on the 
general forms and broad restrictions imposed in the context of a complete (i.e., fully associative) potential 
structure. In particular, in this formulation there are two key ingredients: (1) the thermodynamic 
potential, or Gibbs function, defining the equations of state, and (2) the viscoplastic dissipation function 
for the ensuing flow and evolution laws, of inelastic strain and internal state parameters, respectively. 

In order to progress further, one has to model the specific features of the material concerned (e.g., take 
into consideration experimental observations and also, whenever possible the physics of the phenomena 
involved). This then calls for (1) introducing specialized forms, in terms of the state and internal- 
parameter sets selected, of the potential functions mentioned previously and (2) a careful investigation 
of the impact of these forms on the characteristics of the resulting constitutive model. A detailed study 
involving these two points is the main focus of the present paper. To this end, we shall consider a 
fairly- comprehensive set of equations that are developed here as a subclass of what we, in reference 1, 
termed the decoupled-form, that is, a form encompassing such important features as nonlinear-kinematic 
hardening and both dynamic (strain) and static (thermal) recovery effects. For convenience, we shall 
make use of a nondimensional format of a specific model and construct the plots of the associated state 
space [2,3] to gain further insight into the consequences of the parameters and functional forms chosen 
with regard to the model’s predictive capabilities. 

Although several quantitative results and comparisons are given, we still emphasize that the pro- 
posed equations by no means suggest a specific model intended for immediate utilization. This is 
particularly true with regard to the cyclic response that was modeled (especially in the presence of 
nonlinear-kinematic hardening) where measures such as irreversibility constraints, loading-unloading 
conditions, limit-surface concepts, and so on are often utilized to prevent abnormalities in the predicted 
stress-strain curves upon unloading and reversed loading. We have introduced such internal state dis- 
continuity inequalities in the present model by following approaches currently in common use [e.g. 4-6]. 
However, more work is certainly needed for a valid assessment of such an adhoc treatment, especially 
in the present context of potential-based forms, where the associativeness of the evolutionary equations 
is compromised. Finally, in view of the results obtained here, we firmly believe that the present forms 
provide several important ingredients needed in developing successful models for specific materials. In 
fact, this work is currently underway, and the results for characterization of models of this type will be 
reported subsequently. 

3 Complete Potential Structure 

Here the basic thermodynamic framework put forth by Arnold and Saleeb [1] is summarized. Expressions 
for the Gibb’s thermodynamic and the complementary dissipation potential functions, in terms of a 
number of state and internal variables characterizing the changing internal structure of the material, are 
assumed. For conciseness, the discussion is limited to a case involving small deformations (in which the 
initial state is assumed to be stress free throughout), an initially isotropic material, isothermal conditions 
and the decoupled potential framework [1]. A Cartesian coordinate reference frame and index notation 
are utilized (repeated Roman subscripts imply summation). 

Given the Gibb’s potential in the following form 

$ = $(<Tij,a y ,T, e(j) ( 1 ) 

and assuming a priori that the inelastic strain is an independent parameter (and not an internal state 
variable) for example, 
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( 2 ) 


* = - a ii e{ i + H K) - Z(T) - ^(T - T 0 ) 

an expression for the total strain rate can be obtained by differentiating: 

d,-d$^_ -d*$ A , ./ , f Vrb rx\ 

C °' " dt ( dtTij } " 5cr r , r * + ij +6ii 3 T ( } 

The three terms in equation (3) may then be identified as the reversible, irreversible (inelastic), and 
thermal expansion components of the total strain rate, respectively. Thus, 


where 


and 


p. . — f # _L ff. 4. *T. 


.n _ • 

f,; d<Tijd<T r , <Tr 


■T — c U.T 
€ ij — 3-^ 


(4) 

(5) 

( 6 ) 


but e(j (the inelastic strain rate) is defined separately by using the concept of a dissipation potential 

1 a 7> T). 

Given 

Q = n(<7ij,a 7 ,r) (7) 

then by using the Clausius- Duhem inequality, the flow law can be written as 


. / dn 

€ij ~ dan 


the evolutionary laws for the thermodynamic conjugate internal state variables written as 

an 


Ay — ~~ 


6a* 


and the constitutive equations for the rate of internal state variables as, 

d 7 = [<? 7 *] -1 At 


( 8 ) 

(9) 

( 10 ) 


where 


Qyt = 


-a a H 

da y da 1 


( 11 ) 


Thus equations (8), and (9) represent the flow and evolutionary laws, for an assumed 0 = n(<Tjj , ar 7 , T) 
and equation (10) the internal constitutive rate equations, given a Gibb’s potential $, wherein both 
potentials are directly linked through the internal state variables a 7 . 
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4 Elastic/ Viscoplastic Constitutive Model 

A complete multiaxial statement can be derived by using the above framework, given a specific form 
for both the Gibb’s potential, and the complementary dissipation potential, ft. Form invariance 
(objectivity) of these potentials requires that they depend only on certain invariants of their respective 
tensorial arguments (i.e., an integrity basis [7] ). In the spirit of von Mises and because of the deviatoric 
nature of inelastic deformation, only the quadratic invariant will be considered at this time in specifying 
the dissipation potential. Similarly, only the linear elastic strain energy contribution will be considered 
in specifying the Gibb’s potential, with the internal state groupings being functions of the respective 
quadratic invariants. Finally, although equation (7) indicates that an unlimited number of internal state 
variables can be specified, here our attention will be restricted to a constitutive model with, at most, three 
internal state variables. The first internal state variable, a,j, is taken to be a second order symmetric 
traceless tensor that represents the internal (or back) stress associated with dislocation motion. The 
second and third internal variables (*; and y) are scalar state variables representing the drag and yield 
stresses, respectively, which are associated with dislocation density. 

Consequently, the Gibb’s potential may be written as 

* = -ic rjt( <r r ,<7*, - + tfi(G) + H 2 (V) - Z{T) - - T 0 ) 

and the dissipation potential as 

Q = »J f(F)dF + ^ j g(G)dG + J x{ K )d K +± J z(Y)dY 

where 



(14) 

II 

ps 

(15) 

t = H 

(16) 

K 


Hi(G) = 

(17) 

H 2 (*) = 

(18) 

T 3 

h = 

(19) 

J 2 — -EijSjj 

(20) 

= Sij — a tJ 

(21) 

0 1 c 

Sij — (Tij — 

(22) 

an = atij - 

(23) 


Note that in the preceding expression for the dissipation potential, the stress dependence, both 
external and internal, enters through the scalar functions F and G in the form of effective (E,j) and 


( 12 ) 

(13) 
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internal ( a,*^ ) deviatoric stresses, respectively. Furthermore, the function F acts like a threshold surface, 
because when F < 0, no inelastic strain can occur. Note that the function F was chosen in such a way 
that it would emulate the features of three types (or classes) of viscoplastic theories commonly found 
in the literature. The first class is composed of the models without a yield (or threshold) surface, that 
is, those in which Y = 0 (i.e., $ = 0), which includes models proposed by Miller[8], Walker [9] f Krempl 
et al.[10], and others. The second are those with no evolving drag strength (i.e., k = constant ), yet 
which have a threshold surface, as used by Chaboche[ll], Freed and Walker[12], and others; and the 
third class are those theories with a threshold surface but wherein both the yield and drag stresses are 
proportional, (i.e., Y a k) as put forth by Perzyna[13], Robinson [5, 14] and others. 

By selecting the preceding scalar functions, a general yet complete potential-based model, with asso- 
ciated flow and evolutionary laws that will encompass a wide range of the models found in the literature 
can be constructed. The second invariants , and I 2 , are also scaled for tension. These invariants 
could just as easily have been scaled for shear by replacing the coefficient 3/2 with 1/2, and modifying 
the definition of the magnitude of the inelastic strain rate that follows. Also, by including an additional 
term, /?(?, which is a scaled function of the internal deviatoric stress and drag stress, allows for the 
presence of a strain induced recovery term in the evolution of the associated conjugate variable (*4jbi), 
as discussed later. It is interesting to note the close similarity between the present assumed form of 
the threshold surface F of equation (14) to that used by Lemaitre and Chaboche [15] in the context 
of inviscid plasticity. The significant difference between the two forms is that Lemaitre and Chaboche 
[15] utilized a non- associated format; the form similar to equation (14) was taken to affect only the 
directions of inelastic and internal state whereas a classical form was employed to control yielding 
(i.e., the consistency condition). In the present fully associative case, the same threshold surface F is 
used to affect both yielding,, and the direction of inelastic and internal state. 

Taking the appropriate derivatives of both $ and (i as indicated in equations (2) through (12), one 
obtains the multiaxial isothermal specification particular to the present constitutive model. Here the 
decomposition of the total strain rate is that of equation (4), where the reversible strain rate is, 



efj — CijkiO’ki 

(24) 

and the irreversible (or 

inelastic) strain rate is defined as 




(25) 

where 

11411 =^=^^ 

(26) 

and the internal state variable constitutive laws relating the associated rates of the thermodynamic 
variables and their respective conjugates are 


«ij = [Q 1 tjJfcf] 1 Akl 

(27) 


k = Q2~ l K 

(28) 


Y = Q3 -1 }> 

(29) 

with Ql, Q2, and Q3 being distinctly different operators. 

The evolutionary laws for the thermodynamic conjugate variables are 
(9): 

then obtained from equation 
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( 30 ) 


./ 3 /? .I ./ || 3 Rfi g(G) 

Akl ~ £ *' _ €<i afc,_ ^<c J 


ic — —{Sijffj — aijAij y y } ^g K 


x(k) 


(31) 


And the respective compliance operators from equation (11) are 

Qlijkl = [i aticijj + (33) 

Q1 = M£±il(2p + 3)G’ +I + i(}+ (34) 

Hk* k* 


Q 3 = 


jli 

B K* x ' 


with Iijki = jjt Finally, inversion of the compliance operators results in 


(35) 


[QltjJbl] 1 = 


k 2 H 


3/i(p+ 1 )Gp 


[Ujki - 


3p 


« 2 (2p + 1)G 


akiaij] 


(36) 


Q2" 1 = 


HBk 2 

2 flp(p + l)(2p + 3 )Gp+! + </(« + l)f/*» 


(37) 


Q3- 1 = (38) 

An examination of the evolution of the conjugate of the internal stress (t.e., the internal strain, .4,-j), 
clearly shows that equation (30) possesses, as typically assumed in the literature [16,17], a competitive 
process of a hardening term (which accounts for strengthening mechanisms) and two recovery terms 
(which account for softening mechanisms). The first recovery term evolves with plasticity; it is strain 
induced, and is commonly called dynamic recovery. The second term, interchangeably called static or 
thermal recovery, evolves with time and is thermally induced. 

The evolution of the conjugate of the drag stress (z.e., the internal drag strain, /C), equation (31), 
is driven by three energy sources. The first is the dissipated plastic work (Siji(j) } and the second and 

third, respectively, are the energy changes due to internal structure (dislocation) rearrangement (ciijAij) 
and the yield related dislocation density contribution ( Y y). The static recovery terms ( z(Y ) and x(k)) 
also evolve with time and are thermally induced. Finally, from equation (32) one can see that a scaled 
plastic strain (||cf ||) is driving the evolution of the conjugate of the yield stress (i.e., the internal yield 
strain, y) and is in direct competition with a thermal recovery term z(Y). 


4.1 Specializations 

The preceding model, which consists of three internal state variables, can be specialized into the three 
classes of viscoplastic models discussed previously. Note that in all specializations the form of the 
inelastic strain rate does not change, only its functional dependence on F. Similarly, the applicable 
internal rate constitutive laws (eqs. (27) to (29)) remain unchanged as well. 

Consider the first class of models, those models without a threshold surface ( i.e., no yield stress; 
Y = 0). Here 
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( 39 ) 


F=(£-Mg) 

and the associated evolutionary laws are equations (30) and 

JC — ““ a u -^0 } ” *( K ) (^O) 

which can also be rewritten as 


* = ; { E «4 + + ? ^ 2g } - *w 


(41) 


The second class of models, those containing a threshold surface wherein the drag stress is held 
constant (k = k 0 ), is obtained by assuming 


h±Ph V 2 \ 
*2 K l / 


(42) 


with the associated evolutionary laws given by equations (30) and (32). If we further specialize this model 
to one that has only a single scalar, internal state variable, Y (by eliminating the internal deviatoric 
stress; i.e., 1 2 — 0), equation (42) becomes 



and equation (32) is modified to become 


(43) 


* = < 44 > 

with = | SijSij. Equations (43) and (44) are then similar to other yield strength models found in the 
literature, the only difference being the scaling parameter in the hardening term and the presence of the 
compliance operator Q3 in the internal variable constitutive law of equation (29). 

Finally, the third class of models, those models that possess a threshold surface but which assume 
that both the yield and drag stress are proportional (Y oc «), can be represented by 


F = (S +/?G_1 ) (45) 

with the associated evolutionary laws given by equations (30) and (40). 

Here, the specific form to be investigated will be an internal (or back) stress model of the third class 
(Y = k = k 0 ); that is, 


F = (46) 

with the flow law given by equations (25) and (26) and the internal constitutive law and the evolutionary 
law given by equations (27), (30) and (36). This model, hencforth, will be called the modified model, 
since it is the corrected form of the Robinson model [4,14] as discussed in reference 1, with an additional 
strain-induced recovery term introduced by including the PI 2 term in the threshold function F. The 
primary objectives of this study are to illustrate the effect of the compliance operator (Qlijki) and to 
examine the importance of including a strain-induced (dynamic) recovery term, under both multiaxial 
and uniaxial states of stress, with the limitations imposed by considering only a single tensorial internal 
state variable. This will be accomplished by mapping out the associated state space and the stress-strain 
response. 
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5 Nondimensional Multiaxial Form 

For generality and computational convenience, let us nondimensionalize the specific model to be ex- 
amined by using the following dimensionless variables and parameter groups. Normalizing the applied 
and internal stress components relative to the initial drag strength k 0 of the given material gives the 
dimensionless stresses 


and 


rr . — 11L 
<y X j — 

Ko 


Otij - — 
K 0 


(47) 


(48) 


The time is nondimensionlized relative to the hardening material parameter H , which is defined in units 
per time, by multiplying by the real time t. Thus, 


T H t (49) 

where T is the nondimensional time. As a consequence of equation (49), the rate of change of real and 
nondimensional time are related by 


J_ _ 1 d_ 
dT - Hdt 


(50) 


Henceforth d/dT will be denoted by an open circle appearing over a quantity, that is, (°). Two dimen- 
sionless parameter groups are taken as 


and 


_ R 
n ~ H 

(51) 

B= * 

Hk 0 

(52) 


where ft is a measure of the relative importance of the thermal recovery and hardening terms ( R and 
H ) which appear in the evolutionary equation. Similarly, B may be considered a measure of the relative 
importance of viscous dissipation to plastic dissipation. As pointed out by Robinson [14], B is related to 
the reciprocal of the Bingham number as defined by Prager [18]. 

Given the preceding relations, the nondimensional multiaxial representation of equations (27), (30), 
(36), and (46) can be written as 


and 


o/ 3 || ./ ,1 

if > 0 

(53) 

J 

0= 0 

if Eijdi> < 0 


II o' 

1 1 € ij II - ^ 

!\ iiUi = 2Bf(F)y/7, 

(54) 

0 

« 0 

= [Qlijjti] - 1 Aki 

(55) 


Ju»=2l, -|-^|| || a t , - 3 KBg{*)a kt (56) 
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where 


3e(p + i)S. [;u ‘ , '(2„ + i)c a *'“ iil 

(57) 

and 


F = J 2 + PG - 1 

(58) 

G = h 

(59) 

f 3. . 
h - 2 a i) a i) 

(60) 

J 2 — - E,, 

(61) 

Ejj = Sij &ij 

(62) 

K=(G-G 0 ) Ey[Sijaij] + 6 0 

(63) 

Now further specializing the above expressions by assuming nonlinearity of the inelasticity and ex- 
pressing the thermal recovery dependence by a power law function, we assume 

f(F) = (f)" 

(64) 

and 


g( K) = N m 

(65) 


In the preceding expressions, Orepresents the Macauley bracket operator and Hv[ ] denotes the Heav- 
iside unit function. 

Furthermore, we impose the internal state discontinuity regions suggested by Robinson [5,14]. These 
regions are implied by the resetting of G to an initial value G 0 whenever the product SijOij becomes 
negative in R, equation (63). The creation of these regions is motivated by the essential need to capture 
the cyclic response indicated experimentally within the framework of nonlinear kinematic hardening. We 
emphasize that this is an adhoc treatment in the present context; however, this inconsistency relative to 
the potential structure becomes irrelevant for all of the monotonic- load applications considered in the 
subsequent sections. Hence all conclusions reached remain valid. 


5*1 Uniaxial Simplification and State Space Representation 

If we consider only a uniaxial state of stress, the dimensionless deviatoric applied and internal stress 
tensors become 


and 


Sij 


1 

3 


2d- 0 0 

0 -a 0 
0 0 -<r 



2d 0 0 

0 — d 0 

0 0 -d 
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such that 


F = (<r — d) 2 + Pot 2 - 1 

G = a 2 

and the flow and evolutionary laws, respectively, become 

t =2B{F) n {a -a) (66) 

and 

l - 1 / 1 ff a — \ 1 — (67) 

2fl(p+l)(2p+l)\ (*-o)Jrp (p+l)(2p+l) 

Utilizing the concept of the state space representation [2,3,19], we can examine and illustrate geo- 
metrically the various features of the present unified viscoplastic model, in particular the effect of the 
strain-induced recovery term. The pertinent state spaces are constructed by combining and rewriting 
equations (4), (24), (66), and (67), in standard vector form as follows, and plotting the corresponding Q 
fields or relaxation trajectories. That is, 

S= G(*,&)+V(e) 
or 


£{i- 

where 




«o 

H= — 

2B(p+l)(2p+l) 


(p+l)(2p+l) 

K = (d a — dj)Hv[^da] + a 2 

and e is the nondimensional total uniaxial strain rate. These relaxation trajectories are obtained from 
equations (66) and (68) by imposing (1) a constant nondimensional total strain rate (i.e., e= 0), thereby 
implying that S = and (2) a set of baseline material parameters as follows: 

k = 0.5774 
/i = 1.39xl0" 9 
H = 3.55x10-“ 

R = 3.74xl0 -3 
n = 4 
p = 0.575 
m = 4.365 
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5.1.1 Dynamic and Thermal Recovery Effects 

Figures 1 to 6 illustrate the state space representation associated with equations (66), (68), and the 
base-line material parameters given previously, for three values of and two values of 7 Z . These are, 
respectively; H = 10.535 with /? = 0, 0.5, and 1.0; and 7Z = 0 with /? = 0,0.5 and 1.0. Note that 71 = 0 
can be associated with low homologous temperatures ( T m = < 0.2), whereas It = 10.535 can be 

associated with some elevated temperature {e,g. } T m > 0.35) at which thermally activated mechanisms 
exist . 

Immediately a number of qualitative features of the first quadrant of the state space can be deduced 
from equations (66) and (68). For example, the locus of steady state (or fixed) points in creep, denoted 

o 

by line ab in Figs. 1 to 6, is obtained by imposing the condition S = 0 on equations (68). Mathematically, 
such a condition is obtained when 


<7 = a + 


f . ,y2m+l 'I 

-^ 2 + {* [d _ ( i + /?)d ]} 


l/n-l »/2 


(69) 


and is denoted as the locus of steady state creep, see line ab. The associated steady state creep rate is 

(70) 


0 1_ R d 2m+1 


Q ^ 

Secondly, the locus of maxima in stress (i.e., & = 0) during relaxation (e = 0) occurs when 


0 


or when 

•/ <7 

€ = = 0 
E 

(71) 


3 

“V 

1 

& 

II 

O 

(72) 

This is satisfied when either 

dr — a 

(73) 

(denoted by line de in Figs. 1 

to 6) or 



<j = d + (1 — /Ja 2 ) 1 / 2 

(74) 

(denoted by line ac in Figs. 1 to 6). The points on lines ac and de do not correspond to fixed points in 
the state space, as do those associated with line ab , since the internal stress continues to evolve, that is, 


l = -R d 2("*-p)+i 

(75) 


Note that the inelastic strain rate is zero in this region of the state space as a result of invoking the 
load-unload inequalities in equation (53). Clearly this continued evolution occurs only if the temperature 
is sufficiently high to induce thermal recovery («.e., Tt ^ 0), as in Figs. 1 to 3. For example, when the 
homologous temperature is low (i.e. , T m < 0.2, room temperature for most superalloys), 7Z = 0 and the 
locus of maxima in stress (lines ac and de) will become fixed points (see Figs. 4 to 6). 

Furthermore, comparing the elevated temperature behavior (Figs. 1 to 3) to the low temperature 
behavior (Figs. 4 to 6), we make the following observations. First, with no strain-induced recovery (i.e., 
ft — 0),the model will predict a purely nonlinear hardening response with no secondary creep regime and 
no internal recovery occurring below the line ac, (see Fig. 4). Alternatively, at elevated temperature both 
secondary creep and internal recovery are predicted. Second, when strain-induced recovery is present 
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(:.e., 0 ^ 0) and T m is still low (72 = 0), secondary creep can occur. However ,^no internal recovery will 

take place below the line acde } owing to the inclusion of the inequalities on e . Also, observe that for 
72 = 0 the locus of steady state is linear, whereas for 72 ^ 0 it is a nonlinear function of internal stress 
(a) (see line ab in Figs. 1 to 6). 

Using the evolutionary law in equations (68) and considering a case wherein strain-induced recovery 
is removed (0 = 0), we can show (see Fig. 7) the effect of increasing the thermal recovery term 72 (i.e., 
increasing the temperature) on the relaxation trajectory labeled fa in Fig. 1. Note that 72 has no effect 
on the initial relaxation rate, but rather moves the location of the stationary point (steady state creep) 
up and to the left as indicated by the movement of point g. Furthermore, if 72 = 0, g becomes not 
only a fixed point in creep, but one in relaxation as well. Also for any nonzero 72 value we see that the 
internal stress can, given sufficient time, return to zero, whereas the stress supported by the material 
will return to the initial threshold stress k c . 

Assuming the presence of strain induced recovery mechanisms (i.e., 0 ^ 0, see Fig. 8), we see that 
even when thermal recovery is not present (i.e M 72 = 0) a stationary point in creep will now exist (see 
point g in Fig. 8). Also the line ac and de of fixed points in relaxation become curved and shifted, 
respectively, as compared to Fig. 7. As in figure 7, again we see that including ! 3 and varying 72 does 
not affect the initial relaxation rates (path / to g) until, the point g associated with steady state creep 
is approached and traversed. The initial relaxation rate is, however, greatly affected by changes in the 
nondimensional parameter B as clearly shown in Fig. 9 where B is varied by two orders of magnitude 
while 0 and 72 are held fixed at 0 and 10.535, respectively. Note that because of this change in initial 
hardening, the location of steady state moves along the locus specified by equation (69) since the applied 
and internal stress reached increases with increasing values of B. 

In summary, from the viewpoint of the state space representations the most notable features of the 
present modified model with combined thermal and strain induced recovery are 

1. A finite, elliptical elastic zone that encloses the origin and whose shape is rationalized and suggested 
by Onat in his extensive discussions on state space [3,19] 

2. In creep, the inclusion of strain-induced recovery with an appropriate selection of 0 allows for the 
occurrence of steady state creep outside the elastic zone, even without thermal recovery. 

3. In relaxation, at low temperatures (72 = 0), the present form of the strain induced recovery term 
results in an asymptotic stress value greater than k 0 that is reached within a shorter time interval. 
(Note, Robinson’s original model, with thermal recovery only, will always saturate to a value equal 
to k 0) given sufficient time.) 

5.2 Uniaxial Tensile and Cyclic Response 

Figures 10 and 11 show the predicted tensile response of the present modified model if we assume no 
strain recovery and three different rates of loading (i.e., e= 0.001, 0.01 and 0.1/ hr.) for a low homologous 
temperature (i.e., 72 = 0; see Fig. 10) and a high homologous temperature (i.e., 72 = 10.535; see Fig. 
11). Notice that at a high homologous temperature the tensile response flattens and gives the appearance 
of perfect plasticity, whereas at room temperature the response continues to harden within the range of 
stress investigated. Also, as the rate of loading increases so too does the apparent elastic stress range. 

Introducing strain recovery (i.e., 0 ^ 0) and maintaining 72 = 0 (e.#., room temperature) causes 
the response once again to flatten or saturate to a given level of applied stress (see Fig. 12). This 
clearly illustrates the proposed modified model’s tremendous flexibility in modeling the tensile response 
over a wide temperature range. Note that the ’’knee” of the stress-strain curve is much smoother when 
strain-induced recovery is present (cf. Figs. 11 and 12). Almost identical stress-strain histories are 
predicted, except for the slightly lowered asymptotic stress level, when both strain and thermal recovery 
mechanisms are assumed to be active (see Fig. 13); as compared to the case when only strain-induced 
recovery (Fig. 12) is active. However, comparing the case in which only thermal recovery mechanisms 
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are active (Fig. 11) to that in which there is both thermal and strain-induced recovery (see Fig. 13), we 
see a significant drop in saturated stress levels. 

Finally, Fig. 14 demonstrates the cyclic behavior of the model if thermal recovery and varying 
levels of strain-induced recovery (i.e., /? = 0, 0.5, 1.0) are assumed. Clearly, as the contribution of 
strain recovery is increased, the maximum attainable stress level is reduced. This is consistent with the 
monotonic tensile response predicted earlier, 

5.3 Multiaxial Tension-Torsion Response 

Returning to the multiaxial statement of the proposed modified model given in equations (53) to (63) 
and assuming that the strain-induced recovery term is zero (i.e., /? = 0), one can quickly conclude 
that the only difference between the present modified model and that put forth previously [4,14] is the 
presence of the compliance operator (QUjkl)- The question now becomes whether this operator will have 
a significant impact on the stress-strain response under multiaxial load histories (recall that the uniaxial 
response is unaffected by the presence of this compliance operator). To this end, we will examine the 
effect of this operator under specialized multiaxial load cases, the first one consisting of a proportional 
tension- torsion history and the second, a nonproportional tension-torsion history (see schematic in Fig. 
15). In both cases, stresses will be applied until a specified value is reached and then that stress state 
will be held fixed for a given length of time. 

Considering a proportional load path (i.e., f = X&) one can show both analytically and numerically 
that an identical strain response will be generated for both Robinson’s model [4] (i.e., eq. (36) with 
the second term discarded) and the present modified model (Fig. 16). Figure 16 shows the magnitude 
of inelastic strain (as defined by eqs. (26) and (54)) as a function of nondimensional time T for both 
proportional and nonproportional histories. Note that both Robinson’s model and the current modified 
model give different nonproportional responses, in that the model put forth by Robinson predicts a higher 
accumulation of inelastic strain than does the modified form. This difference is expected, as indicated 
by examination of equation (57). Furthermore, as one might expect, the nonproportional response of the 
two models is initially identical, diverging only when a significant amount of internal deviatoric stress is 
accumulated with respect to the product This point is reached at approximately T= 0.00075 or 

t =2.11 hours. 

At T= 0.001, both <t and f are now constant, and we see that both models predict less accumulation 
of strain than does the proportional case. If, however, the magnitude of the target J 2 is doubled from 
7.2 to 14.4 (cf. Figs. 16 and 17) then Robinson’s model predicts a greater accumulation in inelastic 
strain during loadup then does the proportional case, whereas the modified form still predicts less than 
that of the proportional case. Also, as the load level increases so too does the difference between the 
predicted responses of the two models. For example, at the onset of the constant stress field of J 2 = 7.2 
at T= 0.001 in Fig. 16 and J 2 = 14.4 at T= 0.002 in Fig. 17, the difference in accumulated inelastic 
strain ,||c^||, is 0.0008 and 0.011, respectively. 

To further illustrate the effect of the compliance operator on the evolution of internal state, in Fig. 
18 we show the individual internal deviatoric stress components as a function of nondimensional time 
for Robinson’s model and the present modified form given by equation (55). Note that since R = 0 and 
both models are subjected to the same nonproportional history, the difference in internal state can be 
attributed only to the effect of the compliance operator QUjkl- From Fig. 18 it is clear that the internal 
tensile stress of Robinson’s model evolves while the torsional component remains unaffected , whereas in 
the present modified model the internal torsional stress component decreases until the external torsional 
stress (f) is applied. 

Finally, if strain-induced recovery is included (i.e., /? / 0), the above trends remain; however, the 
rate of accumulation of inelastic strain, in both proportional and nonproportional loading histories, is 
significantly increased (cf. Figs. 16 and 19) because of the material’s internal softening due to the 
mechanical recovery mechanisms. 

If we consider the interaction of (3 and the compliance operator QUjkl relative to the nonlinear hard- 
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ening response, two observations are in order. First, in the uniaxial case either one of these mechanisms 
can easily be fit to reproduce the nonlinear tensile response. On the other hand, for multiaxial load 
applications, a unique feature for nonproportionality or path dependency is gained by the inclusion of 
an internal-state-dependent compliance operator Q. 

6 Conclusions 

A complete multiaxial potential based viscoplastic model has been presented. It contains in general, 
three internal state variables (two scalars and one tensor) and both thermal and strain-induced recovery 
mechanisms. The two scalar internal state variables are associated with dislocation density and are 
defined as the drag and yield stress. The tensorial variable known as the internal (or back) stress is a 
second order, traceless, symmetric tensor and is associated with dislocation motion. By merely modifying 
the definition of the driving threshold function F, this general model was shown to emulate individually 
or simultaneously, three distinct classes of theories found in the literature. These classes are (1) theories 
without a yield surface, (2) theories with a threshold surface yet no evolving drag stress, and (3) theories 
with a threshold surface but with both the drag and yield stress taken to be proportional. 

The remainder of the study was dedicated to the assessment of a single tensorial internal state variable 
(internal stress) model, with the idea of examining the impact of including strain-induced recovery and/ or 
the compliance operator. A parametric study of a non dimensional multiaxial form examined the model 
over a wide range of material properties. The important conclusions of that parametric study are listed 
here: 

• In modeling the stress-strain and creep response of metals in the low homologous temperature 
range, the inclusion of strain-induced recovery provided needed flexibility without adversely affect- 
ing the high homologous response. 

• Inclusion of strain-induced recovery affected the stress level at which steady state creep occurred 
for a given internal stress level, the limit stress under both monotonic and cyclic loading, the 
saturation stress in relaxation, and the elastic domain within the state space. 

• Neither uniaxial nor proportional loading histories were affected by the inclusion of the compli- 
ance operator; however, with multiaxiality the response due to nonproportional loading paths was 
significantly influenced. 

In light of the forgoing, an important task to consider as future work is the development of an appro- 
priate means of maintaining the nonlinear compliance operator under more general cyclic load conditions 
without introducing inconsistencies (e.g., internal discontinuities) into the potential framework. 
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Normalized stress, 



Figure 1.- Calculated {/-field (relaxation lines) 
associated with non dimensional uniaxial sim- 
plification, assuming only thermal recovery (i.e., 
/? = 0; ft = 10.535). 



Figure 3.- Calculated {/-field (relaxation lines) 
associated with the non-dimensional uniaxial 
simplification, assuming both thermal and dy- 
namic recovery (i.e., /? = 1.0; ft = 10.535). 



Figure 2.- Calculated {/-field (relaxation lines) 
associated with nondimen sional uniaxial sim- 
plification, assuming both thermal and dynamic 
recovery (i.e., /3 = 0.5 ;ft = 10.535). 



Normalized internal stress, a 

Figure 4.- Calculated £-field (relaxation lines) 
associated with nondimensional uniaxial sim- 
plification, assuming no recovery mechanisms 

(i.e., 0 = 0; 11 = 0). 


17 




Figure 5.- Calculated {/-field (relaxation lines) 
associated with non dimensional uniaxial sim- 
plification, assuming only dynamic recovery 

(i.e., P = 0.5; Tl = 0). 



Figure 7.- Effect on relaxation trajectory of 
varying the amount of thermal recovery ( i.e., 
H = 0, 1., 10.535, 100.) in the absence of any dy- 
namic recovery (i.e., = 0.). 



Figure 6.- Calculated {/-field (relaxation lines) 
associated with nondimen sional uniaxial sim- 
plification, assuming only dynamic recovery 

(i.e., /? = 1.0 ; 72. = 0). 



Figure 8.- Effect on relaxation trajectory of 
varying the amount of thermal recovery (i.e., 
% = 0, 1., 10.535, 100.) in the presence of dy- 
namic recovery (i.e., /? — 1.0). 
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Figure 9.- Effect on relaxation trajectory of 
varying the amount of viscous-to-plastic dis- 
sipation (t.e., B = 6.781xl0~ 7 , xlO -6 , xlO“ 5 ) in 
the presence of thermal recovery only (/? = 0; 
11 = 10.535). 



0 .2 .4 .6 .8 1.0x10-2 

Normalized total strain, e 

Figure 11.- Effect of strain rate on uniaxial 
tensile response, assuming only thermal re- 
covery (ft = 0 ; = 10.535). 


Normalized stress, 



0 .2 .4 .6 .0 1.0x10-2 

Normalized total strain, e 

Figure 10.- Effect of strain rate on uniaxial 
tensile response, assuming no recovery mech- 
anisms present (/? = 0; 71 = 0). 



0 .2 .4 .6 .8 1.0x10-2 

Normalized total strain, e 

Figure 12.- Effect of strain rate on uniaxial 
tensile response, assuming only dynamic re- 
covery (/? = 1.0; H = 0). 




0 .2 .4 ,6 .0 1.0x10-2 -2.5 -2.0 -1.5 -1.0 -.5 0 .5 1.0 1.5 2.0 2.5x10-3 

Normalized total strain, e Normalized total strain, e 

Figure 13.- feffect of strain rate on uniaxial Figure 14.- Uniaxial cyclic response, assuming 

tensile response, assuming both thermal and thermal recovery along with varying levels of 

dynamic recovery (0 = 1.0; H = 10.535). dynamic recovery (0 = 0, 0.5, and 1.0). 



Figure 15.- Schematic of tension - torsion pro- 
portional and nonproportional load paths, ap- 
plicable for both = 7.2 and 14.4. 


.8x10-2 



0 .5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0x10-3 

NondimensionaJ time, Y 

Figure 16.- Magnitude of inelastic strain ver- 
sus nondime nsional time for both proportional 
and non-proportional loading paths (see fig- 
15) with target J 2 value of 7.2, assuming no re- 
covery mechanisms (0 = 0 ; H = 0; B = 6.781x10’^). 
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Model 
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Figure 17.- Magnitude of inelastic strain ver- 
sus nondimensional time for both proportional 
and nonproportional loading paths (see fig- 
15) with target j 2 value of 14.4, assuming no 
recovery mechanisms (/? = 0; H = 0; B = 6.781x10“^). 


■ 1 Modified, a-j-j 

Robinson’s, 

Modified, a-j 2 

— Robinson's, a *| 2 



Figure 18.- Normalized internal state (both 
normal and shear) versus nondimensional time 
for both proportional and nonproportional load' 
ing paths (see fig. 15) with target J 2 value of 
7.2, assuming no recovery mechanisms (/? = 

0 \K = 0,B= 6.781xl0" 5 ). 
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Nondimensional time, Y 


Figure 19.- Magnitude of inelastic strain ver- 
sus nondimensional time for both proportional 
and nonproportional loading paths (see fig. 
15) with target J 2 value of 7.2, but assuming 
dynamic recovery only (/ 3 = 0.5). 
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